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We construct a two-leg ladder dimer model by using two orbitals instead of one in Su-Schrieffer- 
Heeger (SSH) dimer chain and find out a chiral symmetry in it. In this model, the otherwise-hidden 
additional chiral symmetry allows us to define two chiral massless fermions and show that there 
exist interesting topological states characterized by Z©Z and corresponding zero mode edge states. 

Our complete phase diagram reveals that there is an anomalous topologically nontrivial region in 
addition to the normal one similar to that of SSH model. The anomalous region features that the 
inter-cell hopping constants can be much smaller than the intra-cell ones, being opposite to the 
normal region, and the zero mode edge states do not have well-defined parity each. Finally, we 
suggest that this ladder dimer model can be realized in double-well optical lattices, ladder polymer 
systems, and adatom double chains on semiconductor surfaces. 


The discovery of quantum Hall effect in two- 
dimensional electron gas opens a new era in physics [lE 
@]. In such systems, topology is involved in their elec¬ 
tronic phases and phase transitions. In 2005 Kane et al 
discovered the time-reversal-invariant topological 
phase in graphene with spin-orbit interaction, causing a 
boom in researching topological insulator [6l4l2l|. Since 
topological insulators can be classified according to their 
symmetries fl3| . whether an insulating system is a topo¬ 
logical insulator or not is dependent on the symmetries 
of its crystalline and electronic structures. Fortunately, 
optical lattices engineered with interfering laser beams al¬ 
low us to study specific interesting models which cannot 
be easily found in naturefl4|. Such models can possess 
novel configurations with potentials of single or multiple 
periods. Controls of atoms in the s and p orbitals in op¬ 
tical lattices can exhibit many exotic quantum states [15l— 
[l9| . including topological nontrivial states. For instance, 
there is a topological insulator phase in an s-p-orbital 
ladder reduced from a two-dimensional double-well opti¬ 
cal lattice [20j and a topological semimetal in fermionic 
optical lattice pH- 

Simple dimer chain, also known as Su-Schrieffer- 
Heeger (SSH) model, is a simple but important model 
to show typ ical topological property of one-dimensional 
systems (22|, |23j], although it was originally proposed to 
describe one-dimensional polyacetylene|24|]. The sim¬ 
ple dimer chain belongs to BDI class according to its 
symmetry[l3|, and is characterized by the topological in¬ 
variant C £ Z. 

Here, we construct a two-orbital dimer model, or a lad¬ 
der dimer model, which can be realized by optical lattices 
similar to those with s and p orbitals |15lI20l I2ll] . We find 
that the ladder dimer model can possess higher symme¬ 
try than the dimer chain and thus have richer topological 
phases, beyond BDI class, and can be characterized by 
topological invariant C £ Z © Z. We plot the complete 
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phase diagram and work out its zero-mode edge states. 
Furthermore, we find out the relationship between this 
ladder dimer model and massless chiral fermion. More 
detailed results will be presented in the following. 


Results 


Model and symmetry. We construct a general two- 
leg ladder dimer model by using two orbitals instead of 
one in the simple SSH dimer model. Presented in Fig. 1 
is a schematic of the ladder dimer model. Its Hamiltonian 
can be defined as 

H =^t s (a) B s a jAs + h.c.) + t' s (a) Bs a j+hA s + h.c.) 
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+ tp( a jBp a jA P + h.c.) + t' p (a] Bp a j+1>Ap + h.c.) 

+ ti{ a jBp a jAs + h-C-) + t'i(. a ) Bp aj+i,As + h.c.) 

+ ti{a' jBs a jAp + h.c.) + t' 2 (a i jBs a j+i,Ap + h.c.) 

(1) 

where the fermion operators dip q and a\ Pq (P = A and B , 
and q = s and p x ) are used to define the model. Because 
we are considering s and p x orbitals, with the x axis being 
along the chain, the real-space symmetry of the wave 
functions implies that the nearest-neighbor hybridization 
and hopping can be nonzero and the onsite hybridization 
of the two orbitals needs to be zero. In addition, no 
next-nearest-neighbor hopping is considered to keep high 
symmetry, which is consistent with the SSH dimer model. 
Defining c]=[a] As , a] Bs , a] A? , a] Bp \, we make a Fourier 
transformation under periodic boundary condition, and 
obtain C\ = jC]e ika i = [a{ As , a{ Bs , a{ Ap , al Bp ], 

where L is the number of the unit cells. As a result, 
the Hamiltonian can be simplified into H = E,_ C' k HkCk 
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FIG. 1. A schematic of the ladder dimer model with the 
hopping parameters shown. 


with Hk being expressed as 


/ 0 
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( 2 ) 
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FIG. 2. The reduced energy spectra of E = 
±>/(F + cos fc) 2 + sin 2 fc with £ = 0.5 (upper panel) and £ = 2 
(lower panel), in comparison with the curves described with 
function E = ±(£ + cos k). 


Using two sets of Pauli matrices, = ( t x ,t v ,t z ) in 
the sp orbital space and ~c? = (a x , a y , a z ) in the AB 
subspace, we can express the Hamiltonian in a compact 
form. Because we are considering real parameters, there 
should be some relationships among the parameters [l5l 
[ 20 L [2l| | . If we let the m orbitals be in an order like 
antiferromagnetic order [15J, we can have t\ = t 2 = S 
and t'i = t' 2 = S' for an optical lattice. As a result, the 
Hamiltonian can be simplified as 

H% =1 g [(t + t! cos ka)a x + t' sin kaa y ] 

+ t z ® [(r + t’ cos ka)a x + t’ sin kacr y \ (3) 
+ t x <8> [(<5 + S' cos ka)a x + S' sin kaa y } 

where t = (f s + t p )/ 2, t' = (t' s + t' p )/ 2, r = (t s - t p )/ 2, 
and t' = {t' s — t' p )/ 2. Each of the these parameters can 
be nonzero. Consequently, this model has time reversal 
symmetry (T = K ), particle-hole symmetry (C = I ® 
<j z K ), chiral symmetry (S = I<g>a z ), and space inversion 
symmetry (R = I ® a x ). It belongs to the BDI class and 
is characterized by Z. 

Introducing the unitary transformation defined by U = 
exp(— iaT y /2)®I, we change H% into the following hamil- 
tonian. 

% s k =U^H k U = (t + t' cos ka)I ® a x + t' sin kal ® a y 
+ [(u> + w' cos ka)r z ® a x + w' sin kar z g) a y ] 

(4) 


space, but fortunately, we have only four freedoms in 
two less than those in 'H k . This parameter condition can 
be satisfied by requiring S/S' = t/t' = (t s — t p )/(t' s — t' p ) 
in r H k . It is easy to prove that this hamiltonian com¬ 
mutes with P = t z g) I, [H%,P] = 0. Therefore, we have 
found a hidden chiral symmetry in H k . It should be 
pointed out that the hidden chiral symmetry is guaran¬ 
teed by both the inversion symmetry and the parameter 
condition. With this additional P chiral symmetry, the 
hamiltonian can be viewed as massless fermions and we 
have four additional symmetries: T' = PT , C' = PC, 
S’ = PS, and R' = PR. Then, we can introduce the 
projection operators 

Pl = ^{I + t z )®I,P r = ^(I-t z )®I (5) 

and write the hamiltonian as a block-diagonal form 


H 


S _ 

k 


Hl(k) 0 \ 

0 H^(k) ) 


( 6 ) 


where the H/ = [(t + w) + (f + w') cos ka]a x + (t' + 
w') sin kaa v and H ^ = [(t — w) + (t' — w') cos ka)a x + (t' — 
w') sin kaa y describe the left and right chiral fermions, 
respectively. Each of H/(k) and H^(fc) has two indepen¬ 
dent parameters, and therefore they are independent of 
each other. 

Topological states and phase diagram. Diagonal¬ 
izing the hamiltonians H k and Hr, we obtain the energy 


where w = 5/ sin a and w' = 5'/ sin a, and a (in [ 0 , 7 r]) 
is determined by tana = S/t = S'/t'. Here, the con¬ 
dition 5/t = S'/t' is required in the diagonalization of 
the r space, which reduces one freedom in the parameter 


E 7 = ±yJ Al+Al ( 7 ) 

where 7 = ( L, R), Ai = ( t' ± w') sin ka, A 2 = (t ± w) + 
(t 1 ± w') cos ka, and the minus sign corresponds to the 
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FIG. 3. The ( l,9',ip ') phase diagram for Z| sin (0')|/| sin (<p')| < 
1 for l = r/r' < 1 (a) and 1/Z = r'/r < 1 (b). The true (Z, d,<p) 
diagram is determined by Z| sin (6 + 7r/4)/sin (ip + 7r/4) | < 1 
for the left and Z| sin (6 — 7r/4)/sin (ip — 7t/ 4)| < 1 for the 
right, with 0 < l < +oo. Because of the absolute value 
sign, the complete phase diagram can be plotted in terms 
of 8 G [—7t/2, 7t/ 2] and i p G [0, 7r]. 


left and plus the right. The energy band structure of 
the left part becomes gapless if we have 11 = 1, and 
for the right part the condition is 1| = 1. Each of 
the two parts is a dimer chain. We present in Figure 

2 a reduced form of E 1 , E = ± \J (£ + cos k ) 2 + sin 2 k, 
with £ = = (t + w)/(t' + w') for the left part and 

£ = £r = (t — w )/ (t' — w') for the right one. It is easy to 
see the energy structure inversion in the figure. 

For the hamiltonian H 1 (k) = h^(k) ■ a 

here, the winding number of vector Zi 7 (fc) is 
a good approach to topological characteriza¬ 
tion because there is no a z component [2bl l27j . 
hj(k) = p 7 (fc)(cos 0 7 (fc), sin^ 7 (fc), 0). The winding 


number can be expressed as 

C. = ^ / W* = ^ (8) 


It should be pointed out that our C 7 here is equivalent to 
either zero or +1, and the value -1 does not appear be¬ 
cause we have the same coefficients for sin ka and cos ka 
in il 7 (fc}, which is different from a generalized SSH chain 
model[28j. The total winding number is equivalent to 
C = Cl ® Cr. This definition is well defined because 
it is gauge invariant. Under this definition, the left and 
right parts enter nontrivial topological states respectively 
if the following conditions are satisfied. 


I&-I = 

, t + W 

1 < 1, 

for the left 

1 + 

ee. 

\€r\ = 

1 < 1, 

for the right 

't' — w' 


Generally speaking, the two parts become trivial or non¬ 
trivial independently, and therefore the topological in¬ 
variant belongs to Z © Z. 

The original four parameters, namely t , w , t' , and w ', 
are not convenient to describing the topological phase 
diagram effectively. Instead, we use r, r', 9 , and ip in 
terms of the definition: t = rsinZZ, w = rcosZZ, t' = 
r' sin<p, and w' = r'cosip. Because t' = w' = 0 means 
that the model is trivial, we can assume r' > 0 and r > 0. 
The conditions ([9|) of the nontrivial topological states are 
equivalent to 


. sin (9 + 7t/ 4) sin (9 — 7t/4) 

l\— -7 - L rk \< 1 or ^ h—7 - ITT < 1 

sm (ip + 7T/4) sin [ip — 7t/4) 


( 10 ) 


where l = r/r'. The sign of absolute value guarantees 
that the period is 7r for both 9 and ip. For brevity, we 
show a (1,9' , ip') phase diagram for l\ \ < 1 in Figure 

3. We have two phases with C 7 = 0 and C 7 = 1 for 9'. 
The true phase diagram of the left part can be obtained 
by replacing 9' with 9 + tt/ 4, and that of the right part 
with 9 — 7t/4. The whole phase diagram consists of four 
phases with C = (0 © 0), (0 © 1), (1 © 0), and (1 © 1). It 
can be proved that there is no C = (OffiO) phase for l < 1 
and no C = (1 © 1) phase for l > 1. It should be pointed 
out that the (0 © 1) phase is topologically different from 
the (1 © 0) phase. This is because Hl(k) is indepen¬ 
dent of Hft(k) in H s k . Furthermore, it can be proved 
that there exist a pair of zero mode edge states at the 
interface between the (1 © 0) phase of (£l,£h)=(£i,1/£i) 
and the (Offil) phase of (£l,£r)=(1/£i,£|i), where £i < 1 
is assumed. Therefore, the four phases are topologically 
different from each other. 

It is interesting to compare this phase diagram with 
the conventional dimer model. Here r describes the hop¬ 
ping parameter within the unit cell, and r' between the 
unit cells. The conventional dimer chain is topologically 
nontrivial when and only when the hopping within the 
unit cell is smaller than the hopping between the unit 
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cells. In contrast, our ladder dimer model can be topo¬ 
logically nontrivial even when r is larger than r'. Fur¬ 
thermore, we find that for r < r' , the left and right parts 
can be topologically nontrivial at the same time, but for 
r > r' , they will never be topologically nontrivial simulta¬ 
neously. These happen because for each of the two parts 
the hopping terms are determined by two independent 
parameter freedoms (not independent parameters). 

Zero mode edge states. To explicitly elucidate the 
nontrivial topological states, we study the band struc¬ 
tures and explore edge states under open boundary con¬ 
dition. As unveiled by Delplace et al [26| , under the 
open boundary condition, a finite single chain with L c 
dimers possess two edge states located at the two ends 
when v/v' < 1 — 1 /(L c + 1), where v is the intra-cell 
hopping parameter and v' the inter-cell one. For the 
ladder dimer model, the conditions of edge states are 

l \v+^\ < 1 - CT for the left and l \j^\ < 1 - T7+I 

for the right. In long chain limit of L c approaching to 

oo, these edge-state conditions are consistent with bulk 
conditions for nontrivial topological states. For arbitrary 
finite L c , we define V = l/[ 1 — 1 /{L c + 1)] and obtain the 
conditions for the existence of the zero mode edge states, 


sinfS + tf/l) or smjS 
sin (p + 7t/4) sin (p 


tt/4) 


r/4) 


< 1 , 


( 11 ) 


which takes the same form as the bulk conditions in Eq. 
(10). Therefore, the phase diagram of the zero mode edge 
states can be obtained by substituting l' for l in Eq. (10) 
and Figure 3. This makes a clear correspondence between 
the bulk nontrivial topological states and the zero mode 
edge states. It should be pointed out that there is a 
special region determined by p- < | gin(y±g/ 4 ) l < I near 
the phase boundary, in which there is no zero mode edge 
state but there exist bulk nontrivial topological states. 
The zero mode edge states will appear in pair because 
we have two edges. For l 1 > 1, there can be a pair of 
zero mode edge states in some ( 9 : p ) region and there is 
not any edge state in the remaining region. For l' < 1, 
there are either one or two pairs of zero mode edge states 
for given (9,p) values. The difference between l and l 1 
can be attributed to different finite size effect under the 
open boundary condition from that under the periodic 
boundary condition. 

In order to show the edge states more clearly, we 
present in Figure 4 the energy levels as functions of the 
parameter p with given parameters of L c = 20, r = 1, 
r' = 2, and 9 = Here, because of l = 1/2 and 

l' = 21/40, it can be inferred from Eqs. (10) and (11) 
that at least one of the two parts are topologically non¬ 
trivial, and there are either 2 or 1 pairs of zero mode edge 
states. For the p values near both p = +7r/4 and — 7r/4, 
there are one pair of zero mode edge states with opposite 
parity. In contrast, there are two pairs of zero mode edge 
states near p = 0, with one pair of them belonging to the 
left part and the others to the right part. Each of these 
edge states has well-defined parity. The transition points 
are determined by | sin(<^ ± 7t/4)| > y/3l'/2 = 21^/3/80. 



FIG. 4. The energy levels as functions of the parameter p 
with L c = 20, r = 1, r' = 2, and 9 = 57r/12. At least one 
of the two parts is topologically nontrivial because of l = 1/2 
and l' = 21/40. The inserts (al) and (a2) are the two zero 
mode edge states for p = 7t/4 (— 37t/4 also), corresponding to 
the nontrivial right part. The inserts (cl) and (c2) are the two 
edge states for p = —7r/4 (37t/4 also), corresponding to the 
nontrivial left part. The (bl-b4) shows the four edge states for 
p = 0 (±7t/2 and ±7r also), which corresponds to the situation 
in which both of the two parts are topologically nontrivial. 
The edge states have well-defined parity and appear in pair. 



FIG. 5. The energy levels as functions of the parameter 9 
with L c = 20, r = 2, r' — 1, and p = 5iv/12. In this case 
of l = 2 and l' = 21/10, at most one of the two parts can be 
nontrivial. The inserts (al) and (a2) are the two zero mode 
edge states for 9 = — 7r/4 (37t/4 also), corresponding to the 
nontrivial right part. The inserts (bl) and (b2) are the two 
edge states for p = tv/4 (—37t/4 also), corresponding to the 
nontrivial left part. The edge states appear in pair and have 
no defined parity. 


Presented in Figure 5 are the energy levels as functions 
of the parameter 9 with L c = 20, r = 2, r' = 1, and 
p = y^7 r. In this case, we have l = 2 and l 1 = 21/10, 
and there are at most one pair of zero mode edge states 
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according to Eq. (12). It is clear in Figure 5 that there 
is no edge state near 9 = 0 and equivalent points, and 
there are one pair of zero mode edge states near 9 = 
q=7r/4 and equivalent points. Because being degenerate 
in energy, they have no well-defined parity, in contrast 
to those in the case of V < 1. The transition points are 
determined by | sin(0 ± 7 t/ 4)| < 1/(2/') = 5/21. The two 
edge states can be transformed into each other by the 
space inversion. 


Experimental realization. Now we address how to 
realize the ladder dimer model in (1). Inspired by earlier 
workJTH, [20j . we propose to realize the one-dimensional 
model with a two-orbital optical lattice with the optical 
potential given by 

V(x, y ) = V x i sin 2 kx + V X 2 sin 2 (2 kx + ^) 

l ( 13 ) 

+ V yl sin 2 ky + V y2 sin 2 (2 ky + ^), 


Discussion 

Zak phases. In addition to the winding number, the 
Zak phase can be used to characterize the ladder dimer 
model. The Zak phase for 7 ( L,R ) can be defined as 

$2 = dp{v? q \d q v? q ), (12) 

where the | u^) are the Bloch wave functions, with 7 being 
both L and R. For such one-dimensional models, is 
defined up to 27r due to the gauge transformation. For 
each of L and R, the Zak phase is quantized so that = 
0, 7 r Mod[27r], because of the inversion symmetry. Such 
Zak phase has been measured in the case of simple SSH 
model[29|. Defining Z 1 = 4 >](/ 7 r, we have Z 7 = 0 and 1 
accordingly. It can be proved in terms of definitions ( 8 ) 
and (12) and the gauge invariance that Z 7 is equivalent 
to C 7 . As a result, we can also use the Zak phases to 
characterize the topological phases of the model. 

High chiral symmetry and Z © Z uniqueness. At 
first sight there would be many high-symmetry ladder 
dimer models from the SSH dimer mode, but actually it 
is not the case. First of all, both the on-site hybridiza¬ 
tion and the next-nearest-neighbor hopping needs to be 
forbidden to keep the T, C, S, and R symmetries and the 
hidden chiral symmetry P. This requires that any ver¬ 
tical rung is absent in Fig. 1, in contrast with multi-leg 
extension with vertical rungs [261]. In addition, there is a 
restrictive parameter condition between the inter-orbital 
and the intra-orbital hopping parameters. The hidden 
chiral symmetry is necessary to the additional T', C', S', 
and R' symmetries and the resultant Z ® Z topological 
insulators in the ladder dimer model. The two condi¬ 
tions for nontrivial topological phases in the ladder dimer 
model, < 1 for the left part and < 1 for the right 
in (9), are each similar to that in the SSH dimer model. 
The unitary transformation U allows one more param¬ 
eter freedom in the original hamiltonian H£. Actually, 
one cannot achieve higher chiral symmetry in such two- 
leg ladder dimer model from the SSH model. Therefore, 
in this sense, our ladder dimer model possessing the high 
chiral symmetry (including the hidden chiral symmetry) 
and the Z © Z topological insulator phases is unique up 
to a unitary transformation. 


where we keep V y i, 2 ^ I4i, 2. This optical lattice has 
a double-well structure in the y direction and causes the 
particles to transit alternately between the sub-wells with 
two different tunneling barriers (B\ and B 2 ) in the x 
direction. The tunneling barrier between the two wells 
in the y direction is much larger than B 1 and B 2 in the 
x direction. Consequently, the low-energy physics of the 
two-dimensional system (13) reduces to a two-leg ladder 
model with alternate barriers between sub-wells in the 
x direction. We can control parameters <f> and V y i ?2 to 
tune the well depth of the two legs, and change 141,2 to 
tune the depth of the sub-wells in each leg. Letting the 
s orbital of the upper leg have the same energy as the 
p x orbital of the lower leg, we realize the ladder dimer 
model showed in Fig. 1 and defined in Eq. (1) by filling 
the low-lying levels of the relevant s and p x orbitals with a 
single species of fermions [THUS- Making the p a orbitals 
be in an order like antiferromagnetic order [15J, we can 
obtain t\ = t 2 = 8 and t\ = t' 2 = 8' and simplify the 
Hamiltonian (1) into Eq. (3). In addition, this model 
should be realized in ladder polymer systems or adatom 
double chains on appropriate semiconductor surfaces. 

Conclusion. In summary, we have constructed a 
ladder dimer model with higher symmetry, defined in 
Eq. (3), than the usual SSH dimer chain model. In 
this case, there is a hidden chiral symmetry between the 
two orbitals that allows us to define two chiral massless 
fermions. Furthermore, we show that the ladder dimer 
model can exhibit interesting topological states charac¬ 
terized by Z © Z and have zero mode edge states, assum¬ 
ing a one-dimensional Z ® Z topological insulator. We 
find out complete phase diagram for the bulk topologi¬ 
cal states and zero mode edge states. Importantly, the 
phase diagram reveals that there exist nontrivial topo¬ 
logical states in not only the normal region with r' > r 
but also the anomalous region with r' < r, where r' is 
the inter-cell hopping constant and r the intra-cell one, in 
contrast with the usual dimer chain model which shows 
nontrivial topological states only in the r' > r region. 
We also find that the zero mode edge states have well- 
defined parity in the normal region, but do not in the 
anomalous region. Finally, we suggest that this ladder 
dimer model can be realized in double-well optical lat¬ 
tices, ladder polymer systems, and adatom double chains 
on semiconductor surfaces. 
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